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abstract 

The masses of supermassive black holes correlate almost perfectly with the velocity dispersions of their host 
bulges, M. oc cr", where a = 4.8 ±0.5. The relation is much tighter than the relation between M, and bulge 
luminosity, with a scatter no larger than expected on the basis of measurement error alone. Black hole masses 
estimated by Magorrian et al. (1998) lie systematically above the M. - a relation defined by more accurate mass 
estimates, some by as much as two orders of magnitude. The tightness of the M m - a relation implies a strong link 
between black hole formation and the properties of the stellar bulge. 



1. INTRODUCTION 

After decades of indirect and circumstantial evidence, the 
motion of gas and stars on parsec scales has provided irrefutable 
dynamical evidence for the presence of 10 7 - 10 9 M black 
holes (BHs) in about a dozen elliptical and a handful of spi- 
ral galaxies (Kormendy & Richstone 1995). While efforts to 
build a larger, statistically significant sample continue, we have 
now moved from debating the existence of supermassive BHs, 
to asking what regulates their formation and evolution, and how 
their presence influences, and is influenced by, their host galax- 
ies. 

In an early review based on eight detections, Kormendy 
& Richstone (1995) found that BH masses M, scale linearly 
with the absolute blue luminosity of the host bulge or ellipti- 
cal galaxy. This correlation was later strengthened by Magor- 
rian et al. (1998) using a larger (~ 30) sample of galaxies to 
which simple stellar dynamical models were applied. At the 
same time, it has been noted (e.g. Jaffe 1999) that the M, —B^ 
relation suffers from observational biases and exhibits a large 
scatter which is not accounted for by the uncertainties in the 
individual measurements. 

By understanding how the properties of BHs relate to those 
of their host galaxies, we can hope to learn about the formation 
and evolution of both. In this Letter, the connection between 
BH masses and the stellar velocity dispersion of the host galaxy 
is investigated for the first time. We find a remarkably tight 
correlation with negligible intrinsic scatter when using galax- 
ies with well-determined BH masses (roughly speaking, those 
galaxies in which the observations have resolved the sphere of 
gravitational influence of the BH). Our results suggest that the 
stellar velocity dispersion may be the fundamental parameter 
regulating the evolution of supermassive BHs in galaxies. 

2. database 



All secure BH mass estimates available to date (see §3), to- 
gether with a compilation of properties of the host galaxies, are 
given in Table 1 . Revised Hubble and T-type (from the Third 
Reference Catalogue, RC3, de Vaucouleurs et al. 1991) are 
found in columns 2 and 3, while column 4 lists distances to the 
host galaxy. With a few exceptions detailed in the footnotes, 
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all distances are from surface brightness fluctuation (SBF) data 
(Tonry et al. 2000) calibrated as in Ferrarese et al. (2000). 

Total apparent magnitudes rmj, uncorrected for Galactic ab- 
sorption, are from the RC3 for all elliptical galaxies (T-type = 
-4 or smaller), and from de Vaucouleurs & Pence (1978) for 
the Milky Way. For the lenticular and spiral galaxies (T-type = 
—3 and larger), trig for the bulge is derived using the empirical 
correlation between T-type and the ratio between bulge and to- 
tal luminosity (Simien & de Vaucouleurs 1986), and is deemed 
to be accurate within 0.5 mag. Finally, all magnitudes are cor- 
rected for Galactic extinction using the DIRBE/IRAS maps of 
Schlegel, Finkbeiner & Davis (1998) and an extinction law fol- 
lowing Cardelli, Clayton & Mathis (1989), and converted to 
absolute magnitudes (col. 5) given the distances in column 4. 

The methods used in deriving the BH masses, and references 
to the original papers are listed in the last column of Table 1 . 
Because the masses depend linearly on the assumed distance to 
the host galaxies, the values in column 6 have been corrected 
to adhere to our homogeneous set of distances. This correc- 
tion is random in nature, and negligible with the exception of 
IC 1459, which is twice as distant as assumed in the original 
paper. Uncertainties in the host galaxies' distances have been 
incorporated in the errors in the BH masses. 

Elliptical galaxies and bulges of spirals have radial velocity 
gradients, hence a measure of the velocity dispersion a will de- 
pend on the distance to the galaxy, the size of the aperture used, 
and the location of the aperture with respect to the galaxy core 
(e.g. Davies et al. 1987). For this work, we have chosen the 
same definition of a used for studies of the fundamental plane 
of elliptical galaxies, namely the central velocity dispersion, 
typically measured in an aperture a few arcseconds in diameter 
(Davies et al. 1987, and references listed in the footnotes of 
Table 1). Our choice will be justified in §3. To bring all val- 
ues of a to a common system, we have adopted the prescription 
of Jorgensen, Franx & Kjaergaard (1995) and transformed all 
velocity dispersions to the equivalent of an aperture of radius 
r e /8, where r e is the galaxy (or bulge) effective radius. The ap- 
plied corrections are very small (rarely exceeding 5%) and are 
deemed accurate to within 1% (Jorgensen, Franx & Kjaergaard 
1995). Raw and corrected a are listed in columns 7 and 8 of 
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Table 1. Database of secure black hole mass estimates and properties of the host galaxies (Sample A) 



Galaxy 


Rev. Hubble 


T-type 


Distance^ ^ 




BH Mass 


<T< 2 ) 




Method and 


Name 


Type 




(Mpc) 


(mag) 


(10 8 M@) 


(km s- 1 ) 


(km s- 1 ) 


References^ 3 ) 


MW 


Sbl-II 




8.0±0.9(') 


-19.13±0.50 


0.0295±0.0035 


100±20 


100±20 


PM; GeOO 


11459 


E3 


-5.0±0.3 


30.3±4.0 


-21.50±0.32 


4.6±2.8 


322±41 


312±41 


G; VKOO 


N221 


cE2 


-6.0±0.3 


Q.8±Q.l 


-15.76±0.18 


0.039±0.009 


80±10 


76±10 


S; vdM98 


N3115 


scr 


-3.0±0.3 


9.8±0.6 


-19.74±0.52 


9.2±3.0 


291±38 


278±36 


S; EDB99 


N3379 


El 


-5.0±0.3 


10.8±0.7 


-20.03±0.14 


1.35±0.73 


210±27 


201±26 


S; GOO 


N4258 


SAB(s)bc 


4.0±0.3 


7.2±0.3 


-18.26±0.51 


0.3901±0.034 


146±19 


138±18 


M; M95 


N4261 


E2 


-5.0±0.3 


33.0±3.2 


-21.26±0.22 


5 4+ 1 ' 2 

° -1.2 


306±40 


290±38 


G; FFJ96 


N4342 


scr 


-3.0±0.5 


16.7±1.0 


-17.24±0.52 


3.3+j-; 


255±33 


261±34 


S; CB99 


N4374 


El 


-5.0±0.3 


18.7±1.2 


-21.44±0.15 


17 +12 
1 '-6.7 


304±39 


286±37 


G; B98 


N4486 


EOpec 


-4.0±0.3 


16.7±1.0 


-21.61±0.14 


35.7±10.2 


370±48 


345±45 


G; M97 


N6251 


E 


-5.0±0.8 


104±10 


-21.94±Q.28 


5.9±2.0 


293±38 


297±39 


G; FF99 


N7052 


E 


-5.0±0.7 


66.1±6.4 


-21.33±0.38 


■> 7+2.6 


270±35 


261±34 


G; MB98 



1 For lack of independent determinations, distances to NGC 6251 and NGC 7052 are derived as v/Hq, where the systemic 
velocities v are from the CfA redshift survey (Huchra et al. 1998) and H = 71 ± 7 km s _1 Mpc -1 (Mould et al. 2000). The 
distance to NGC 4258 is geometrically determined from the proper motion of its nuclear water masers (Herrnstein et al. 1999). 
NGC 4342 has been assumed to be at the same SBF distance as the nearby NGC 4472. The distance to the galactic center is in 
kpc (from Genzel et al. 2000). 

2 A11 velocity dispersions are from Davies et al. (1987), except for the Milky Way (Kent 1992), NGC 4258 (Terlevich et al. 
1990), NGC 6251 (Heckman et al. 1985) and NGC 7052 (Wagner et al. 1988). 

3 Codes for the original papers are given in the references. Codes for the methods used in estimating the masses are: G = gas 
kinematics from HST optical spectra; S — stellar kinematics from HST optical spectra, using axisymmetric dynamical models 
with 3-integral distribution functions; M = kinematics of water maser clumps, derived from VLBA data; PM = proper motion 
measurements of the SgA star cluster. 



Table 1 . For the Milky Way, we adopt the mean velocity dis- 
persion compiled by Kent (1992) within r e /%. 

3. ANALYSIS 

We did not make any attempt to homogenize the error es- 
timates on the BH masses. Except for adding in quadrature 
the (small) uncertainty in the galaxies' distances, the errors on 
M, listed in Table 1 are those quoted by the respective au- 
thors. However, the real uncertainties are often much larger. 
For instance, Magorrian et al. (1998) derived BH masses based 
on fitting a simple class of dynamical models to ground-based 
kinematical data. In almost all of these galaxies, the data can 
equally well be fit by a more general class of model with no 
BH at all; thus the Magorrian et al. mass estimates might con- 
servatively be interpreted as upper limits (e.g. van der Marel 
1997). The same is true for the majority of ground-based, 
stellar-kinematical BH detections. 

In view of this fact, we list in Table 1 only the galaxies which 
we deem to have reliable BH mass estimates. The proper mo- 
tion studies of the Sagittarius A star cluster and the dynamics 
of the water maser disk in NGC 4258 lead to the most robust 
determinations of M.. Close seconds are estimates in 10 addi- 
tional galaxies, based on data from high-resolution HST obser- 
vations, either absorption-line stellar spectra or observations of 
the motion of nuclear dust/gas disks. These 12 galaxies com- 
prise our "Sample A." Additional galaxies with less secure BH 
mass detections are Arp 102B (M, obtained from fitting accre- 
tion disk models to variable optical emission lines, Newman et 
al. 1997), and all galaxies for which M, is estimated based on 
stellar kinematics obtained from the ground (e.g. the Magorrian 
et al. 1998 sample). These galaxies define our "Sample B": the 
data are tabulated in Merritt & Ferrarese (2000). 



We then searched for linear correlations between logM. and 
both Bj and logtr c . We used the bivariate linear regression rou- 
tine of Akritas & Bershady (1996), which accommodates intrin- 
sic scatter as well as measurement errors in both variables. M. 
was taken as the dependent variable. The results of the regres- 
sion fits, applied to each sample of N galaxies, are summarized 
in Table 2 and Figure 1 . 

The correlation between M, and bulge magnitude (Fig. 
la,c) is poor, both for Sample A and Sample B. Although the 
best linear fit to the data has a slope close to the value of -0.4 
expected if M, is simply proportional to the bulge mass, it is 
apparent from the figure, and from the reduced x 2 r of the fit (Ta- 
ble 2), that even by restricting the sample to the galaxies with 
the most accurately determined BH masses, the intrinsic scat- 
ter in the M % -B { \ relation remains significantly larger than the 
reported errors. No sub-sample of galaxies, selected either by 
Hubble Type or method used in deriving M., defines a tight lin- 
ear relation between M, and Bj . This implies that differences 
in the mass-to-light ratio between Hubble types, or systematic 
biases affecting any particular method, are unlikely to account 
for the large scatter. 

Figures lb and Id show the dependence of M, on the cen- 
tral stellar velocity dispersion a c of the host bulge or elliptical 
galaxy. The correlation is remarkable: Sample A, which shows 
a large scatter in the M, -Bj plots, now defines a linear relation 
with negligible intrinsic scatter. The best-fit linear relation is 

logAf. =4.80(±0.54)logcr e -2.9(±1.3) (1) 

with M, in units of M Q and a c in km s" 1 . The slope of the 
relation remains unaltered, albeit with a larger uncertainty, if 
the two galaxies at the low-velocity-dispersion end of the dis- 
tribution (the Milky Way and M32) are excluded from the fit. 
The reduced \ 2 of the fit (Table 2) is only 0.8, consistent with a 
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scatter that derives entirely from measurement errors. The first 
incarnation of Eq. (1) was suggested by Merritt (2000) (the 
"Faber-Jackson law for black holes"). 

The galaxies in Sample B define a much weaker correlation 
between M, and a c (Fig. Id). Furthermore, the BH masses in 
this sample lie systematically above the mean line defined by 
Sample A, some by factors of ~ 10 2 . Two factors distinguish 
the two samples: the reliability of the M. estimates; and the 
method used to derive the BH masses. About 1/2 of the mass 
determinations in Sample A are based on gas motions while al- 
most all of the Sample B masses are derived from stellar kine- 
matics. We see no evidence for a systematic difference between 
the two types of mass determination; for instance, NGC 4342 
and NGC 7052 have identical a c and M., even though the de- 
termination of M, in NGC 4342 is based on stellar kinematics 
and in NGC 7052 on rotation of a gas disk. In the case of IC 
1459, for which M, predicted by Eq. 1 is 2.5 times larger than 
measured, Verdoes-Kleijn et al. (2000) suggest that the true 
BH mass could be a factor 3-4 greater than their best estimate 
due to non-circular motions of the gas. It seems likely that the 
different correlations defined by the two samples result largely 
from errors in the determination of M, for the galaxies in Sam- 
ple B. 

Our choice of aperture-corrected, central velocity dispersions 
is convenient but not unique. We note first that correcting a for 
the effect of aperture size does not introduce a bias in either the 
slope or the intercept (see Table 1). However, the need for 
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FIG. 1 . — (a): BH mass versus absolute blue luminosity of the host elliptical 
galaxy or bulge for our most reliable Sample A. The solid line is the best linear 
fit (Table 2). Circles and triangles represent mass measurements from stellar 
and dust/gas disk kinematics respectively. The squares are the Milky Way (M. 
determined from stellar proper motions) and NGC 4258 (M, based on water 
maser kinematics), the only two spiral galaxies in the sample, (b) Again for 
Sample A, BH mass versus the central velocity dispersion of the host elliptical 
galaxy or bulge, corrected for the effect of varying aperture size as described 
in §2. Symbols are as in panel (a), (c): Same as panel (a) but for Sample B. 
Circles are elliptical galaxies, squares are spiral galaxies. The solid line is the 
same least-squares fit shown in panel (a); the dashed line is the fit to Sample 
B. All BH mass estimates in this sample are based on stellar kinematics, (d): 
Same as panel (b) but for Sample B. Symbols are as in panel (c). 



Table 2. Results of the Linear Regression Fits: Y = aX + (3 



X,Y Variables' 1 ) 


Sample' 2 ' 


N 


a 


P 




log(<Jc), log(M.) 


A 


12 


4.80±0.54 


-2.9±1.3 


0.79 


log(<Jc), log(M.) 


B 


29 


5.81±0.43 


-4.6±1.0 


2.3 


log(<r), log(M.) 


A 


12 


4.81±0.48 


-3.0±1.1 


0.61 


log(«, ma ). log(M.) 


A 


10 


4.61±0.79 


-2.3±1.9 


8.0 


B°, log(M.) 


A 


12 


-0.36±0.09 


1.2±1.9 


23 


B°, log(M.) 


B 


30 


-0.48±0.10 


-0.8±2.0 


25 



1 Units are Mq for Mm, km s 1 for a and a c , and magnitudes for B^, 

2 See §3 for a definition of the samples 



aperture corrections could be avoided by using a measurement 
of the rms velocity at some fiducial distance from the center. 
Figure 2 plots M, vs. the rms stellar velocity v rms at r e /4, with 

v rms = \J {a 2 + v 2 / sin 2 0r f /4- Here a and v r are the measured 
stellar velocity dispersion and mean line-of-sight velocity re- 
spectively. A complication with this approach is the typically 
poorly-constrained value of the inclination angle i between the 
rotation axis and the line of sight. Estimates for ; are avail- 
able only for NGC 3115 (Emsellem, Dejonghe & Bacon 1999) 
and NGC 4342 (Scorza & van den Bosch 1998). Neglecting or 
wrongly estimating sin ; will increase the scatter in the relation 
and bias the slope too low, by moving faint, rapidly-rotating 
galaxies to the left in the M. - v rms plane. Nevertheless, for our 
sample A, linear regression fits (Table 2) show that the slopes of 
M, vs bulge velocity are coincident whether a c or v rms is used. 

An interesting question is whether the tight correlation be- 
tween M m and er c might simply reflect the influence of the BH 
on the stellar kinematics of the nucleus. The coincidence of the 
slopes obtained when v rms is substituted for a c is the most con- 
vincing evidence that this is not the case, since v rms is measured 
well beyond the radius at which the BH could have a measur- 
able effect. In addition, most of the measurements of a listed in 
Table 1 were carried out using apertures much larger than the 
expected radius of gravitational influence of the BH. We stress 
that - even if the correlation between M, and a c were due in 
part to the gravitational influence of the BH on the motion of 
stars in the nucleus - this would not vitiate the usefulness of 
the relation as a predictor of M m . Figure lb suggests that M, 
can be predicted with an accuracy of <~ 30% or better from a 
single, low-resolution observation of a galaxy's velocity disper- 
sion. This is a remarkable result. 

4. DISCUSSION 

We have found a nearly perfect correlation between the 
masses of nuclear BHs and the velocity dispersions of their host 
bulges, M m oc a a , a = 4.8 ± 0.5. Here we examine some of the 
implications of this correlation. 

The Magorrian et al. (1998) mass estimates fall systemati- 
cally above the tight correlation defined by our Sample A (Fig. 
Id), some by as much as two orders of magnitude. The discrep- 
ancy is a strong function of distance to the galaxy, particularly 
at the high-mass end: nine of the Magorrian et al. galaxies have 
BH masses that are larger than the largest BH mass in our Sam- 
ple A (3.6 x 1O 9 M , in NGC 4486), and six of these are more 
distant than 50 Mpc. A number of authors (van der Marel 1997; 
Ho 1998) have suggested on other grounds that the Magorrian 
et al. mass estimates may be systematically high. If our Eq. 
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(1) correctly predicts M # , the gravitational radius of influence 
of the BHs in most of these galaxies would be far too small to 
have been resolved from the ground. For example, Eq. (1) pre- 
dicts M. ~ 2.8 x 10 8 M Q for NGC 4874, a full two orders of 
magnitudes smaller than the Magorrian et al. estimate; the im- 
plied radius of influence is ~ 24 pc ~ 0"05. In support of this 
idea, we note that the best-fitting M, found by Magorrian et al. 
in five of their 36 galaxies was negative, while an additional 
three galaxies - altogether, 1/4 of their sample - were consis- 
tent with M. < 0. In view of this, we suggest that correlation 
studies based on the Magorrian et al. masses (e.g. Merrifield, 
Forbes & Terlevich 2000) be interpreted with caution. 




60 80 100 200 400 

(7 or v rms (km s~' ) 



FIG. 2. — BH mass versus the central velocity dispersion a c of the host el- 
liptical galaxy or bulge (solid circles) or the rms velocity v rms measured at 1/4 
of the effective radius (open circles). Crosses represent lower limits in Vrms- 
The solid and dashed lines are the best linear fits using u c (as in Figure lb) and 
v rms respectively. 



In passing, we caution against the indiscriminate extrapola- 
tion of Eq. (1) much below the range plotted in Figure 1 (for 
example to the range appropriate to dwarf elliptical galaxies 
or globular clusters), as the formation mechanism of BHs with 
masses smaller than ~ 10 s M Q might differ from that of more 
massive systems (Haehnelt, Natarajan & Rees 1998). 

Why should BH masses be so tightly correlated with bulge 
velocity dispersions? One possibility is a fundamental connec- 
tion between M, and bulge mass, with a a good predictor of 
bulge mass - a better predictor, for instance, than B®. This 
explanation is superficially plausible, since the masses of early- 
type galaxies scale with their luminosities as M~ L 5 / 4 (Faber 
et al. 1987) and L - a 4 , hence M ~ a 5 . The M.-a relation 
of Figure lb would therefore imply a rough proportionality be- 
tween BH mass and bulge mass, i.e. that a universal fraction 
of the baryonic mass was converted into BHs. However early- 
type galaxies appear to be two-parameter systems (Djorgovski 
& Davis 1987) and it is not clear that a alone should be a good 
predictor of galaxy mass. 

Another possibility is that a measures the depth of the poten- 
tial well in which the BH formed. A number of authors (Silk & 
Rees 1998; Haehnelt, Natarajan & Rees 1998) have suggested 
that quasar outflows might limit BH masses by inhibiting ac- 
cretion of gas. Equating the energy liberated in one dynamical 
time of the bulge to the gravitational binding energy, and as- 
suming accretion at the Eddington rate, gives a maximum BH 
mass that scales as a 5 (Silk & Rees 1998), again consistent with 
the observed relation. This dependence could be maintained in 
the face of mergers only if BHs continued to grow by gas ac- 
cretion during all stages of the merger hierarchy (Kauffmann & 
Haehnelt 2000). 

We thank Avi Loeb for an illuminating discussion which 
served as a key motivation for this paper, and W. Dehnen, C. 
Joseph, J. Magorrian and M. Sarzi for comments and sugges- 
tions which improved the manuscript. LF acknowledges grant 
NASA NAG5-8693, and DM acknowledges grants NSF AST 
96-17088 and NASA NAG5-6037. This research has made use 
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